arXiv:1508.07264v2 [physics.plasm-ph] 23 Nov 2015 


Direct and secondary nuclear excitation with x-ray free-electron lasers 
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The direct and secondary nuclear excitation produced by an x-ray free electron laser when interacting with 
a solid-state nuclear target is investigated theoretically. When driven at the resonance energy, the x-ray free 
electron laser can produce direct photoexcitation. However, the dominant process in that interaction is the 
photoelectric effect producing a cold and very dense plasma in which also secondary processes such as nuclear 
excitation by electron capture may occur. We develop a realistic theoretical model to quantify the temporal 
dynamics of the plasma and the magnitude of the secondary excitation therein. Numerical results show that 
depending on the nuclear transition energy and the temperature and charge states reached in the plasma, 
secondary nuclear excitation by electron capture may dominate the direct photoexcitation by several orders of 
magnitude, as it is the case for the 4.8 keV transition from the isomeric state of ®^Mo, or it can be negligible, 
as it is the case for the 14.4 keV Mossbauer transition in ®^Fe. These findings are most relevant for future 
nuclear quantum optics experiments at x-ray free electron laser facilities. 
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I. INTRODUCTION 

The recently operational X-ray Free Electron Laser 
(XFEL) facilitiesii^ have considerably enhanced the focus 
on x-ray physics, with particular emphasis on x-ray quan¬ 
tum optics^ and the extension of concepts used in atomic 
physics towards the high-frequency radiation domain. A 
peculiar circumstance is that x-rays are no longer reso¬ 
nant to valence-electron transitions in atoms. Instead, 
the corresponding resonant systems are either inner shell 
electron transitions in (highly) charged iona^“— , or tran¬ 
sitions in atomic nuclei^i^. Nuclear quantum optics fea¬ 
turing the interaction of x-ray light with Mossbauer nu¬ 
clei in the few keV transition energy range has gained 
considerable momentum, both theoreticallji^rJ^ and 
experimentallySiiSr^i. Interestingly, these experiments 
operate with less than one resonant x-ray photon per 
pulse on average^^ since they make use of synchrotron 
sources with brilliance values about eight orders of mag¬ 
nitude lower than the ones at the XEEL. This raises the 
question what qualitative and quantitative differences are 
to be expected in nuclear quantum optics experiments to 
be performed in the near future at XFEL sources? Con¬ 
sidering the fundamental laser-matter interaction pro¬ 
cesses involved, this is also a new situation compared 
to photonuclear studies with petawatt optical lasers2^“— 
like, for instance, the currently built LEEX^^ in Osaka 
or ELI-NP— in Bucharest. 

While so far in all x-ray-nucleus experiments the elec¬ 
tronic response only acted as background, the increase 
of the electric field strength leads to drastic changes in 
the interaction between photons and electrons which may 
additionally influence the nuclear excitation. Eor photon 
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energies of approx. 10 keV and middle-range atomic num¬ 
ber Z materials, the photoeffect dominates the electronic 
processes. Due to the unique interaction between high 
intensity x-ray pulses and matter— new matter states 
like cold, high-density plasmas can be reached^i^. In 
such environments secondary nuclear processes from the 
coupling to the atomic shell are rendered possible by the 
presence of free electrons and atomic vacancies. For in¬ 
stance, in the resonant process of nuclear excitation by 
electron capture (NEEC)^, a free electron is captured 
into a bound state simultaneously transferring the ex¬ 
cess energy to the nucleus, as illustrated schematically 
in Fig. [H While the related process nuclear excitation 
by electron transition (NEET) has been observed, e.g., 
in experiments with an experimental evidence 

of NEEC remained elusive so far, mainly due to strong 
radiative background processes. 

Secondary nuclear processes in the XEEL-produced 
plasma open new channels of nuclear excitation. De¬ 
pending on the goal of the nuclear quantum optics ex¬ 
periment, this might be of advantage, for instance if the 
maximization of nuclear excitation is desired, or of dis¬ 
advantage, if coherence-based enhancement of nonlinear 
interaction between x-rays and nuclei is envisaged. In 
the latter case the interactions in the plasma would lead 
to strong decoherence rates. In a recent worb^ we have 
shown theoretically that the secondary NEEC occurring 
in the plasma environment created by the XEEL pulse 
shining on a bulk nuclear solid-state target can exceed 
by orders of magnitude the direct resonant photoexci¬ 
tation. The concrete example studied was the case of 
photoexcitation starting from the isomeric, long-lived nu¬ 
clear excited state of ®^Mo at approx. 2.5 MeV energy. 
The XEEL-induced 4.85 keV excitation may open in this 
case the possibility to retrieve the energy stored in the 
metastable state in a process known as isomer triggering 
or isomer depletion^^— . In this first comparison between 
direct and secondary nuclear excitation processes, the 
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temporal evolution of the plasma states was neglected, 
i.e., NEEC was assumed to take place in the same plasma 
conditions for a period of 100 ps. In reality, the plasma 
undergoes Coulomb explosion during this time and we 
may expect that only a reduced fraction of the assumed 
100 ps is effectively providing the advantageous parame¬ 
ters for secondary nuclear excitation in the plasma. This 
critical approximation in Ref<^ calls therefore for further 
investigation. 

The purpose of the present work is twofold. Eirst, we 
would like to provide a more reliable theoretical predic¬ 
tion on the secondary nuclear excitation by taking into 
account the time-dependent dynamics of the plasma. We 
do so by parameterizing the plasma expansion by means 
of a hydrodynamical mode l^^d^ in the quasi-neutral limit 
and including atomic processes in the plasma with the 
help of the population kinetics model implemented in 
the ELYCHK code^. Employing this improved model, 
we provide here a more accurate comparison of the direct 
and secondary nuclear excitation in the isomer triggering 
process investigated in Ref.—. The second motivation is 
to provide reliable theoretical predictions for direct and 
secondary nuclear excitation in the case of the 14.4 keV 
Mossbauer transition in ^"^Ee. That nucleus has been 
the “working horse” of the nuclear forward scattering 
community at synchrotron source a^^d^ and it is also the 
present candidate for x-ray quantum optics using nuclear 
transitions in thin film x-ray planar cavitie a^d'^4^ . Many 
of the present studies at synchrotron sources have ex¬ 
tensions envisaged with the XFEL as source of stronger, 
nonlinear nuclear excitation. Eor these studies it is vital 
to know whether also additional plasma effects may play 
an important role for the nuclear excitation or for the 
sought-for coherence effects. 

Our results show that, although for the case of iso¬ 
mer triggering in ®^Mo, secondary nuclear excitation via 
NEEC remains dominant by several orders of magnitude, 
this is not at all the case for ^^Ee, where the secondary ex¬ 
citation can be safely neglected. We identify here criteria 
related to the nuclear transition energy, the atomic struc¬ 
ture and plasma conditions that can be used to identify 
whether for a particular nuclear species, the secondary 
processes can be dominant compared to the direct pho¬ 
toexcitation channel. This knowledge is then applied to 
the present nuclear transition candidates starting from 
stable or metastable ground states and energies in the 
operation range of the XFEL. These results are most rel¬ 
evant for the layout of future nuclear quantum optics 
experiments at the XFEL. 

The paper is structured as follows. In Section El the 
theoretical formalisms used to describe photexcitation 
and NEEC in the plasma environment are briefly re¬ 
viewed. The hydrodynamical expansion model used for 
the plasma is introduced here. Section uni presents our 
numerical results on direct and secondary nuclear exci¬ 
tation. A detailed analysis of the plasma condition is 
included. This section concludes with a comparison be¬ 
tween the two excitation channels and a general set of 



FIG. 1. (color online) Schematic illustration of nuclear exci¬ 
tation by electron capture^. A free electron (right-hand side) 
recombines into a highly charged ion with the simultaneous 
excitation of the nucleus (simplified level scheme on the left- 
hand side). In this example, the ®®Mo nucleus is initially in a 
metastable state. The ®®™Mo excitation is induced by NEEC 
into the L shell (right-hand side) with subsequent decay to 
the nuclear ground state (long blue solid and dashed arrows 
in the left panel). The nuclear levels are labeled with their 
total angular momentum, parity and energy (in keV). Repro¬ 
duced with permission from Phys. Rev. Lett. 112, 082501 
(2014). Copyright 2014 American Physical Society. 

qualitative criteria to identify the dominant nuclear ex¬ 
citation channel. The main conclusions of the paper are 
summarized in the last section. 

II. THEORY 

The theoretical formalisms used to address the direct 
photoexcitation and the secondary NEEC process in the 
XFEL-produced cold plasma are briefly reviewed in the 
following. The plasma dynamics is described by means 
of a hydrodynamical model. The effects of the latter on 
the secondary excitation channel are quantified. Atomic 
units h = me = e = 47reo = 1 are used throughout 
Secs. Ill Al and IIIBI 

A. Direct light-nucleus interaction 

The natural and most convenient way to describe the 
direct interaction of the laser electromagnetic field with 
the two-level nuclear system is via the well-established 
density matrix formalism. The detailed description of 
this treatment has been presented elsewher o'^^d^ and we 
restrict ourselves here to only few relevant expressions. 
The time evolution of the density matrix p is determined 
by the master equation^ 

= [Ho + Hi, p] + Cp , (1) 

where the so-called Lindblad operator C includes deco- 
herent relaxation processes of the system like sponta- 
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neous decay or dephasing of the laser field. The system 
Hamiltonian can be separated into the unperturbed part 
Hq and the interaction part Hi. The dynamics of the 
system is then governed by the interaction Hamiltonian 
Hi given by 

Hi = J d^rjn{r) ■ A{r,t) , ( 2 ) 

where c stands for the speed of light and jn represents the 
nuclear current density. The vector potential of the elec¬ 
tromagnetic field denoted by A is described classically. 
Making use of the multipole expansion and applying the 
long-wavelength approximation it can be shown that the 
interaction Hamiltonian Hi is directly proportional to the 
electric and magnetic multipole moments Gllm and Mlm 
which are defined by 

^LM = J (frr^YLM{0,ip)pn{r) , 

J - Mr) ■ (3) 

Here, L denotes the total photon angular momentum and 
M its projection, represents the nuclear charge density 
and stand for the vector spherical harmonics Yjf 
with J = L^. 

The calculation of the nuclear interaction ma¬ 
trix elements then reduces to the evaluation of 
{IeMe\^LM\IgMg) and {IeMe\MLM\IgMg), where the 
ground (excited) nuclear state wave functions are la¬ 
beled by the total nuclear spin Ig{Ie) and its projection 
Mg(Me). Utilizing the Wigner-Eckart theorerrii^ the ma¬ 
trix elements can be related to the reduced transition 
probabilities B which are defined as^ 

B{£iM)LJg^ Q = ^^\{I,\\&LiML)\\Ig)\^ . (4) 

Here, £ stands for electric and Ai for magnetic parity 
transitions, respectively. Since often experimental val¬ 
ues are available for the reduced transition probabilities 
B, one can avoid the limitations of theoretical nuclear 
models. For an electric (magnetic) transition of fixed 
multipolarity L, the interaction matrix element reads^ 

{IeM,\Hi\IgMg) 



X ^^^y^^,, v/24TT^g(g(Af)L,4^ Je) , (5) 

where k denotes the wave number, ojk the angular fre¬ 
quency and Ek the electric field amplitude of the laser, 
respectively, and C(ji j 2 j 3 ',rni m 2 m 3 ) stands for the 
relevant Clebsch-Gordan coefficient. 

The natural nuclear transition width T 0 that enters the 
Lindblad operator in Eq. is usually composed of two 
parts: (i) the radiative decay where the nuclear excitation 


energy is released in form of a photon, and (ii) internal 
conversion (IC) which accelerates a bound electron into 
the continuum due to the nuclear deexcitation. Typi¬ 
cally, To is in the range of 10“^ to 10“^° eV for low-lying 
nuclear transitions which are accessible with the XFEL. 
The large discrepancy between the nuclear width Tq and 
the XFEL bandwidths (currently in the order of several 
eV) is one of the main limiting factors for the direct laser- 
nucleus interaction. However, if solid-state targets are 
considered, there can be a third decay channel, the col¬ 
lective or coherent decay, which may strongly enhance 
the nuclear transition width and hence the laser-induced 
nuclear excitation. This collective decay relies on the 
fact that in solid-state targets composed of Mossbauer 
nuclei, photons can be absorbed and re-emitted without 
recoil. Due to the recoilless scattering it is impossible 
to distinguish which nucleus has been excited such that 
a collective excitation also known as nuclear exciton is 
created^. The formation of the exciton requires that af¬ 
ter excitation, each nucleus decays to the initial ground 
state. Otherwise, processes like decay by IC, nuclear re¬ 
coil or spin-flip which provide (be it even only in princi¬ 
ple) “which-nucleus” information prevent the formation 
of the collective effect. The collective excitation predom¬ 
inantly decays by emitting a photon into the forward di¬ 
rection with an increased decay rate^. By considering 
a very short incoming laser pulse, approximated as a <5 
function in time, the time dependence of the scattered 
intensity immediately after nuclear excitation can be ap¬ 
proximated as^^— 

lit) oc , (6) 

where ^ = janUQL stands for the dimensionless thick¬ 
ness parameter with rig denoting the number density of 
nuclei in a sample of length L and ctr the resonant nuclear 
cross section. The faster decay rate expression (.^ -|- l)ro 
is only valid in relation to the formation of the nuclear 
exciton. For instance, while the excitonic state is typi¬ 
cally populated in nuclear forward scattering experiments 
with ®^Fe, this cannot be the case for ®^Mo. The reason 
for the latter is that due to the isomer triggering nuclear 
level scheme (see Fig. [1}, the final nuclear state after the 
decay F does not coincide with the initial state IS. This 
allows at least in principle a way to discern which nuclei 
have been excited, and prohibits collective effects. 

B. NEEC: microscopic description 

A detailed presentation of the NEEC theoretical for¬ 
malism can be found in Refsi^i^. For completion, here 
we reproduce some relevant expressions that are used 
in our later numerical analysis. The initial state for 
NEEC describes the system composed of the nucleus 
in its ground or a metastable excited state, the bound 
electrons indexed by the quantum numbers aoi the total 
electronic angular momentum quantum number Jg and 
its projection /ig, and additionally a free electron with 
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asymptotic momentum p and free electron spin projec¬ 
tion TTli^, 


I'tpi) = \IiMi,{aoJofJ.o,P'ms)) ■ (7) 

In the process of NEEC, the free electron is captured by 
simultaneously exciting the nucleus to some state |/dMd). 
The electronic capture level is denoted by the quantum 
numbers ad, Jd and /id- With that, the intermediate 
state can be written as 


IV'd) = |7dMd, ttdJdA^d) ■ (8) 


For the description of the electron-nucleus interaction 
we adopt the Coulomb gauge expression 






Pn(rn) 
kn - Tel 


(9) 


where rn,e stand for the nuclear and electronic spatial 
coordinates, respectively. Important to note is that the 
attractive Coulomb interaction in Eq. (0 can only me¬ 
diate transitions of electric multipolarity. The mag¬ 
netic electron-nucleus interaction is instead transferred 
via the exchange of a virtual photon between the nuclear 
and electronic currents^. The corresponding interaction 
Hamiltonian is obtained by performing a perturba¬ 
tion expansion of the system’s transition operator and 
making use of the Feshbach projection method^^. Ac¬ 
cording to Refj^ the magnetic Hamiltonian can be writ¬ 
ten as 

i7„, = --a./dVn^^4^ ’ (10) 


where a is the vector composed of the three Dirac ma¬ 
trices ax, ay and az- 

For the calculation of the NEEC transition rate we as¬ 
sume that the projection quantum numbers of the initial 
and final states are not resolved in the experiment. This 
is accomplished by averaging over M,, /ig and TOs and 
summing over and ^d- Furthermore, we integrate 
over the solid angle Up of the incoming electrons. For 
the general case of both electric and/or magnetic transi¬ 
tions, the NEEC transition rate reads^^*^ 


Txi->-d 
^ neec 


27r 

2(2/i + l)(2Jo + l) 


E E 

Mirria/io MdMd 


X |(V'd|77en +77m|'0i)l^/Oi , 


drip 


( 11 ) 


where pi stands for the initial density of free electron 
states. The final result for the NEEC rate defined above 
for capture into ions with closed shells (Jg = 0), a certain 
multipolarity L and a transition of electric parity then 
reads 


yiSL) 
^ neec 


+ ^ /d)(2Jd + 1) 


E<E c[j,Lr,-,:-r 


( 12 ) 


where j is the total angular momentum of the continuum 
electron, and n is the Dirac angular momentum quantum 
number related to j via j = |k| — The index Kg de¬ 
notes the Dirac angular momentum quantum number of 
the bound electron after the capture. In the case of a 
magnetic transition with multipolarity L we obtain 


^ ^d)(2Jd + 1) ^ 


L2(2L-t !)■ 

X (2j -|- l)(Kd + k)^ 


Jd j L 

1 _i 0 

2 2 


R 


(M) 


(13) 


In the expressions above, Rn stands for the nuclear radius 
and the occurring radial integrals Rf\^ ^ and Ri^^ ^ 
require a numerical approach to be solved. Their defini¬ 
tions can be found in Refs, ^^’^^ . 

The NEEC cross section is then given by 

27r2 

<-:ei{E) = —Ketc Td(E - Ed) (14) 

pZ 

where p is the absolute value of the initial free elec¬ 
tron momentum p. The function Lg represents the well- 
known normalized Lorentz profile occurring in resonant 
systems defined as 


with Ed and Td being the energy and the natural width 
of the resonant state, respectively. For arbitrary capture 
levels ipd the width Eg is composed of an electronic and 
a nuclear part, Eg = E^ -|- E™'”/ 


C. NEEC reaction rates in the plasma environment 

In a plasma, free electrons with different kinetic ener¬ 
gies are available. NEEC on the other hand is a reso¬ 
nant process where a free electron is captured by an ion 
with the simultaneous excitation of the nucleus. The res¬ 
onance bandwidth is determined by the Lorentz profile 
m- Since the kinetic energy of free electrons in a plasma 
is distributed over a wide range, many resonant NEEC 
channels may exist. In the following we will shortly de¬ 
scribe how such a situation can be handled theoretically 
in terms of reaction rates. 

As introduced in Sec. Ill Bl the initial and intermediate 
states are given by Eq. ([7l) and Eq. dS]), respectively. In 
order to restrict the number of possible initial configura¬ 
tions, for a lower-limit estimate, we consider in the fol¬ 
lowing only NEEC into ions which are in their electronic 
ground state. In this case, the initial electronic configura¬ 
tion ag is uniquely identified by the charge state number 
q before electron capture. In the isolated resonance ap¬ 
proximation, the total NEEC reaction rate in the plasma 
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can be written as a summation over all charge states q 
and all capture channels a^, 

A„eec = ^^P,Are^c , ( 16 ) 

q ad 

where the partial NEEC rate into the capture level Od of 
an ion in the charge state q is given by 

A^e^ = J dEal^tiE)ME) ■ (17) 

The single-resonance cross sections are defined in 
Eq. (fldll . The dependence on q is hidden here in the 
NEEC resonance energy Ed entering in the Lorentz 
profile. The factor Pq occurring in Eq. (USD denotes the 
probability that an ion of charge state q is present in 
the plasma. The electron flux rpg in the plasma can be 
written as the product of the density of states g{E), the 
Eermi-Dirac distribution fpD^E, Tg) for a certain electron 
temperature Tg and the velocity v{E), 

ME) = g{E)fME, Tg)v{E) . (18) 

The temperature dependence of 4>e is only included in the 
Eermi-Dirac statistics /fd- The density of states as well 
as the velocity are determined by taking the relativistic 
dispersion relation for the free electrons. The chemical 
potential /ie of the electrons occurring in /fd is fixed by 
adopting the following normalization 

j dEg{E)fMEM) = ng . (19) 


D. Plasma expansion 

In the scenario under investigation, the plasma- 
formation occurs on the time scale of the XFEL pulse 
duration (^ 100 fs) while the plasma expansion time 
is in the range of picoseconds. Accordingly, we neglect 
the plasma expansion during its formation in the laser- 
target interaction. In order to ascertain the effect of the 
plasma expansion on atomic processes after the inter¬ 
action of the laser pulse with the target, we model the 
expansion of the target plasma by a quasi-neutral ex¬ 
pansion of spherical clusters as studied in the context of 
the intense optical laser pulses interaction with spheri¬ 
cal cluster a^^i^^i^^~ — . We follow the analysis of Ref4^ 
to describe the plasma expansion. The target plasma 
is assumed to maintain a uniform (but decreasing) den¬ 
sity throughout the plasma sphere during the expansion 
while the electron temperature decreases with the adia¬ 
batic expansion of the plasma, 

^rigV dTg = -Pg dV , (22) 

where rig is the number density of free electrons, and V 
is the volume of the plasma with the radius R, i.e., V = 
47ri?^/3. Here and in the following, we use the Lorentz- 
Heaviside natural units /i = c = fee = cq = /xo = 1 for 
the plasma modeling part. The pressure of free electrons 
Pg is given by the ideal gas law, Pg = rigTg. Therefore, 
the time-dependent electron temperature Tg{t) and the 
plasma radius R{t) satisfy the relation 

Te = Te,0 , (23) 


Thereby, rig represents the number density of free elec¬ 
trons. By using the definition of the NEEC cross section 
and assuming that the free electron momentum and the 
NEEC interaction matrix elements are constant over the 
width of the Lorentz profile Ld(A — Ed), Eq. (1171) can be 
simplified to 


27r2 

\q,ad _ yi->d^res 
^neec o neec 

pZ 


(Ed) 


where the resonant electron flux is defined by 


( 20 ) 


where Te,o is the initial electron temperature and i?o the 
initial plasma radius. During the plasma expansion, the 
electrons lose their thermal energy to the ions resulting 
into the electron and ion kinetic energies 




= -rie 


_d 

dt 



(24a) 



(24b) 


®r(Ed) = J dE Ld{E - Ed)ME) . 


with mi being the ion mass. The equation of plasma 
(21) expansion then reads^ 


The net NEEC rate Aneec provided by Eq. (fTO is 
strongly dependent on the available charge states and free 
electron energies which are both dictated by the plasma 
conditions. In turn, the latter will not be constant over 
time as the plasma undergoes expansion. We proceed to 
formulate a model that takes the spatial expansion of the 
plasma into account and provides the temporal dynam¬ 
ics of the plasma parameters required to calculate the net 
NEEC rates in the plasma environment. 


d'^R rig Tg o Rl Tg^ 

"''IF = =*y-R5- = 


(25) 


where Z\ = ngjn-i is the ratio of the electron density to 
the ion density. In the quasi-neutral limit, E, is therefore 
the average charge state of the ions. Integrating once for 
a fixed Zi = Eg yields 


dR 

dt 


2 



( 26 ) 
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where Vs = {STefiZo/mi)^^^ is the ion sound velocity 
which, in the limit i? —>■ oo, is the characteristic speed for 
the plasma expansion^i^. It may be noted that Peano 
et analyzed the expansion of spherical nanoplas¬ 

mas with the Vlasov-Poisson equations, the particle-in¬ 
cell (PIC) method, and the ergodic model^. Compar¬ 
isons between their results^^*^ and those via the hydro- 
dynamic expansion [Eg. (1251) ] (atomic processes not in¬ 
cluded) for plasmas with the dimensionless parameter 
To = 3Ayi?g = 7.2 X 10-3 or Tq = 7.2 x lO’^, where Ad 
denotes the Debye length for the electrons, show that the 
hydrodynamic expansion model can adequately describe 
the expansion of spherical cluster targets heated by an 
intense optical laser pulse. In addition, ID and 2D PIC 
simulations using the EPOCH code^ for plasmas (atomic 
processes not included) with the parameters under con¬ 
sideration in this article have also been performed. The 
expansion time obtained by the hydrodynamic expansion 
model is in good agreement with the PIC simulation re¬ 
sults (not shown here). 

The quasi-neutral hydrodynamic plasma expansion de¬ 
scribed above requires the plasma density and the charge 
state as input parameters. One can exploit the different 
time scales involved in the plasma formation and subse¬ 
quent plasma expansion, and as a first step calculate the 
average charge state by using the electron temperature 
from the plasma expansion model. One can then use this 
charge state to estimate the plasma expansion and con¬ 
sequently the electron temperature at a later time, thus 
establishing a feedback loop between the effective atomic 
processes (represented by the average charge state) and 
the plasma expansion. This is discussed in Sec. IIIIEI 


III. RESULTS AND DISCUSSION 

We investigate and quantify the direct and secondary 
nuclear excitation in a normal incidence setup for two 
species of nuclei, (i) ®3]\/[o isomer triggering where the 
isomeric state IS is depopulated via a 4.85 keV transition 
to an above lying triggering level T, as depicted in Fig.[T] 
and (ii) the resonant driving of the 14.4 keV transition 
in 3'^Fe from the ground state G to the first excited level 
E. For both nuclear targets, we assume a thickness of I 
fim, which for the case of ^^Fe corresponds to a resonant 
thickness of ^ = 3.9. The laser focal spot is assumed to 
be 10 /rm^, i.e., focal radius of approx. 1.8 /rm. 

In the case of ^'^Fe the initial state is the stable ground 
state of the isotope, such that bulk iron samples enriched 
with the 3'^Fe isotope (natural abundance 2.2%) can be 
fabricated. The ®3]V[o solid-state sample should on the 
other hand contain nuclei in the isomeric state at 2.5 
MeV excitation energy. The isomers can be produced 
in 4fNb(p,n)®3™Mo reactions^, directly embedded into 
1 fim thick solid-state niobium foils^. Based on the pro¬ 
duction reaction cross section values^ (^ 30 mb), we 
estimate that a isomer density of 10^® cm-^ can 

be achieved in the solid-state 1 /rm-thick Nb foils^ us¬ 


ing standard proton beams like the LINAC at 
(Considering possible transport difficulties between iso¬ 
mer production and final experimental sites, it would be 
of great advantage to have an in situ production of the 
isomers directly at XFEL facilities, in conjunction with 
optical-laser-driven proton acceleration.) The majority 
of atoms in the isomeric sample belong to the niobium 
species {Z = 41) and for the plasma parameter estimates, 
their interaction with the XFEL photons will play the 
dominant role. However, due to the very close atomic 
number values of molybdenum and niobium, the photo¬ 
electric response of the two atomic species is very similar. 

In the following we address in more detail the XFEL 
parameter choice and present our photoexcitation and 
NEEC numerical results for the two considered nuclear 
transitions. 


A. Light sources 

The advent and commissioning of the XFEL promise 
significant progress in the field of nuclear quantum op¬ 
tics. It offers in particular the possibility to investigate 
the direct laser-nucleus interaction with coherent, highly 
brilliant x-ray pulses. Characteristic for this new kind of 
radiation sources are the small wavelengths in the x-ray 
regime, high power, high brilliance and coherent light 
pulses. The pulse duration lies usually in the range of 
hundreds of femtoseconds. 

Currently, there are two operating XFEL facilities 
worldwide, the LCLSi at SLAC in Stanford and the 
SACLA^ in Japan. The LCLS provides photons with an 
energy of approximatively 10 keV with an average spec¬ 
tral brightness of up to 2.7 x 10^^ photons/(s mrad^ mm^ 
0.1%bandwidth). The SACLA facility in Japan delivered 
the so-far highest photon energy of 19.5 keV. In addition, 
several XFEL machines are in construction like the Eu¬ 
ropean XFEL at DESY in Hamburg, the SwissFEL at 
the Paul Scherrer Institute in Switzerland, and MaRIE 
at the Los Alamos National Laboratory in the United 
States. The European XFEL, for instance, is expected to 
achieve photon energies up to 24.8 keV. The XFEL gen¬ 
erated higher harmonics may even provide photon pulses 
with energies above 25 keV. 

In addition to the brightness of the XFEL radiation, 
the coherence properties are in particular prominent in 
comparison to broad-band synchrotron radiation. How¬ 
ever, while the XFEL pulses are fully spatially coherent, 
their temporal coherence is poor because of random fluc¬ 
tuations in the electron charge density at the start-up of 
the SASE process. Until now, there are two ideas how to 
tackle the problem of the poor temporal coherence. The 
first idea is to load the undulator with an already seeded 
light pulse in order to reduce shot-to-shot fluctuations at 
the start-up (seeded XFEL)22i“— . The second idea is to 
construct an x-ray cavity which directs the light pulse 
several times through the undulator (XFEL oscillator, 
XFELO)^^. 
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Parameter 

LCLS 

SACLA 

Eur. XFEL 

XFELO 

Tmax (eV) 

10300 

19600 

24800 

25000 

BW 

2x10“® 

2.2x10"® 

8x10“'* 

1.6x10“* 

Tpuise (fs) 

100 

100 

100 

1000 

Tcoh (fs) 

2 

a 

0.2 

1000 

Ppeak (W) 

1.5-4x10“ 

10“ 

2x10*° 

4.1x10° 

T W \ 

-^peak 

3.9x10“ 

9.8x10“ 

2.0x10*'* 

4.0x10*° 

/rep (Hz) 

30 

10 

2.7x10* 

10° 


® In our calculations we assumed 10% of the pulse duration, i.e., 
10 fs. 


TABLE I. The approximate maximal achievable photon en¬ 
ergy i?inax, bandwidth BW, pulse duration Tpuise, coherence 
time Tcoh, peak power Ppeak, peak intensity /peak and pulse 
repetition frequency /rep for the four considered XFEL facil¬ 
ities: LCLS^-^, SACLA42i2i, Eur. XFEL^ and XFELO^^. 
A focal spot of 10 is assumed. 


In this work we consider the parameters of the cur¬ 
rently running XFELs, the LCLS and the SACLA, and 
of the European XFEL still under construction at DESY. 
Additionally, we will also provide results for the XFELO 
in order to point out the importance of temporal coherent 
radiation pulses in the direct laser-nucleus interaction. 
The corresponding laser parameters are summarized in 
Table m Here, we assume a moderate laser focusing on a 
spot of 10 /rm^ (a focal length of 7 nm has for instance 
already been achieved in Refi^). 

A general hindrance in the laser-induced nuclear exci¬ 
tation is that only a small fraction of the laser photons 
really fulfills the nuclear resonance condition due to the 
usually small nuclear transition widths. This has two 
major consequences. First, since only a small number 
of nuclei is excited per pulse, the repetition rate of the 
XFEL facility is a key ingredient for an effective driving of 
nuclear transitions. The repetition rates of the LCLS and 
the SACLA have with 30 Hz and 10 Hz, respectively, the 
same order of magnitude. The future European XFEL is 
expected to provide light pulses with a repetition rate of 
2.7 X 10'’’ Hz and the XFELO even with a frequency of 
10® Hz, which is 3 and respectively 4-5 orders of magni¬ 
tude higher than the corresponding values of the already 
operational XFELs. 

The second consequence is that we need to introduce 
an effective laser intensity which accounts for the mis¬ 
match between photon and nuclear transition energy^, 
namely 

. (27) 

^ laser 

Thereby, Fnuci denotes the nuclear transition width and 
Fiaser the bandwidth of the laser pulse. The effective 
field amplitude Aes required in Eq. ([S]) in order to obtain 
numerical results is proportional to yjles. 


B. Direct photoexcitation 

The photoexcitation rates are calculated following the 
formalism in Ref4i using reduced transition probabil¬ 
ity values from Refi^. The limited laser coherence time 
is accounted for by introducing a corresponding deco¬ 
herence rate in the Lindblad operator. The numerical 
solution of the master equation for the density matrix 
is then carried out with the procedures implemented in 
Mathematical. We considered realistic parameters for 
the coherent high-frequency light sources shown in Ta¬ 
ble HI Results for the total population of the excited level 
Pgg®®’' after a single laser pulse are presented in Table HIl 
These populations are obtained by summing the elements 
Pee{Me) = (A Afe|p|/eMe) evaluated at t = Tpuise over 
all possible projection quantum numbers Mg. For both 
93mMo and ®^Fe cases, the XFELO delivers the highest 
excitation rates per pulse among the considered laser fa¬ 
cilities. 

Comparing the effective laser intensities displayed in 
Table HD it can be seen that the XFELO provides with 
6.7x10^2 W/cm^ for ^^'"Mo and 4.0x10“ W/cm^ for 
®^Fe the highest value of les- The effective intensity of 
the LCLS, for instance, is 3 orders of magnitude smaller. 
Thus, due to its narrow bandwidth, the XFELO is able 
to provide in average 1000 times more resonant photons 
per pulse than the LCLS, SACLA or the Eur. XFEL. 

Another, possibly less obvious reason for the outstand¬ 
ing excitation properties of the XFELO are its coherence 
properties. As already mentioned in the previous section, 
all currently operating XFELs lack a good temporal co¬ 
herence (indicated by Tcoh in Table HD because of the 
random fluctuations in the initial electron charge den¬ 
sity. The poor temporal coherence is one of the main 
limiting factors for the nuclear excitation. Considering 
for instance totally coherent x-ray pulses (Tcoh = oo) in¬ 
creases the nuclear excitation between 4 and 6 orders of 
magnitude. 

Experimentally, the nuclear excitation is best acces¬ 
sible by measuring the number of photons of a specific 
transition emitted in the decay of the excited nuclear 
state E. In the case of isomer triggering, for instance, the 
triggering level T first decays to an intermediate state F 
which then subsequently falls down via a decay cascade 
to the nuclear ground state G. In this cascade a charac¬ 
teristic photon of 1 MeV is emitted which can be used 
as signature for the isomer triggering. The rate of 7 -ray 
signal photons at the detector is directly proportional to 

Pee I 

fiaser ^ ^ (38) 

^ 1 + aic 

where Nq = ngAfocT represents the number of nuclei 
present in the focal spot Afoc and aic stands for the IC 
coefficient of the transition producing the 1 MeV signal 
photon. Furthermore, denotes the branching ra¬ 

tio, i.e, the probability of a nucleus in T to not fall back 
to the state IS. For ®®Mo, this probability that the trig- 











®3Mo 



°"Fe 


les (W/cm®) 

laser 

Pee 

(1/s) 

7eff (W/cm®) 

laser 

Pee 

(1/s) 

LCLS^ 

5.2x10® 

1.9x10"®° 

5.6x10"^^ 

3.1x10® 

9.1x10"’° 

1.9x10"® 

SACLA 

1.2x10® 

1.7x10"®° 

1.7xl0"^'‘ 

7.0x10’’ 

9.5x10"’° 

6.6x10"® 

Eur. XFEL 

6.5x10® 

2.4x10"®^ 

6.5x10"^® 

3.9x10® 

1.2x10"’° 

2.2 

XFELO 

6.7x10^® 

4.5x10"^"' 

4.6x10"° 

4.0x10“ 

2.2x10"® 

1.5x10® 


In the case of ®^Fe, we assumed the x-ray pulse resonant to the 14.4 keV transition, although i?max = 10.3 keV for the LCLS. 

TABLE II. Excited state occupation number and signal photon rate for ®®Mo and ®^Fe. Both and have 

been calculated by considering only the direct photoexcitation. The laser parameters shown in Table U are employed here. 


gering level decays to the ground state via the emission of 
the I MeV signal photon can be approximated by one. In 
the case of ®^Fe, the resonantly scattered photons via the 
transition from the excited state E to the ground state 
G can be used as detection signal. In this case, a simi¬ 
lar expression as Eq. (1281) can be used, with a branching 
ratio value representing the probability that the excited 
nuclear state E decays radiatively. Taking the collective 
channel into account, this probability is approx. 82% for 
a sample thickness of 1 p,m = 3.9). 

Results for are also shown in Table |lTl Since the 

nuclear excitation per pulse is typically very small, the 
pulse repetition frequency /rep of the laser plays a cru¬ 
cial role in order to have detectable signal rates. The 
Eur. XFEL, for instance, is expected to provide two or¬ 
ders of magnitude larger rates than the LCLS although 
the excited state population per pulse is smaller, sim¬ 
ply due to the high repetition frequency of 27000 pulses 
per second. The XFELO may produce 4.6 x 10“^ and 
1.5 X 10® signal photons per second for ^^'"Mo and ®^Fe, 
respectively. The large difference between ^^'"Mo and 
®^Fe comes mainly from the magnitude of the interaction 
matrix elements and the number of nuclei present in the 
sample, which are orders of magnitude apart. This dif¬ 
ference also compensates that the effective intensity for 
the resonant photoexcitation of ®^Fe is about one order 
of magnitude smaller than for 


C. Secondary NEEC in a stationary plasma 

The net NEEC rate Aneec provided by Eq. mi) is 
strongly dependent on the available charge states and free 
electron energies which are in turn both dictated by the 
plasma conditions. The charge state distribution (CSD) 
can be calculated by applying the collisional-radiative 
model implemented in FLYCHK—. In that model the 
CSD is completely determined by fixing the electron tem¬ 
perature Te and the ion density rij present in the plasma. 
In the following we investigate the role of the plasma con¬ 
ditions, considered to be stationary here, on Aneec for the 
numerical case of ^^'"Mo isomer triggering. 

The microscopic transition rates Fneec given in Eq. m 
are computed numerically following the formalism de- 


capture 

^res 

qlS^F 

‘^neec 

orbital 

[l/m®/s/eV] 

[b eV] 

3^3/2 

1.10 X 10®’ 

3.04 X 10"® 

3^5/2 

1.09 X 10®’ 

4.28 X 10"® 

4^3/2 

2.77 X 10®° 

1.05 X 10"® 

4^5/2 

2.76 X 10®° 

1.50 X 10"® 

5^3/2 

1.59 X 10®° 

5.10 X 10"® 

5^5/2 

1.59 X 10®° 

7.28 X 10"® 


TABLE III. NEEC case study for a 350 eV plasma with ions 
initially in the charge state q = 24. The resonant electron 
flux {Te = 350 eV) and the NEEC resonance strength values 
are presented for captures into d ^/2 and d ^/2 orbitals. 


veloped in Ref4^. The electronic wave functions and 
binding energies are calculated by the relativistic multi- 
configurational Dirac Eock (MODE) method imple¬ 
mented in the computer code GRASP932^. The prob¬ 
abilities Pq are determined by the FLYCHK-CSD results 
and the electron flux (j>e is calculated following the ex¬ 
pressions in Sec. Ill Cl Finally, the convolution integral 
over the Lorentz profile and the electron flux appearing 
in Eq. dig is solved by approximating the resonance pro¬ 
files by Dirac delta distributions centered at E^- 

In Eig.Ha) NEEC cross sections are displayed for sev¬ 
eral capture levels together with the electron flux Eor 
the latter we have exemplarily taken the temperatures 
Te = 350 eV and Te = 500 eV. The figure shows that the 
NEEC cross section decreases by going to more loosely 
bound capture levels. Due to the reduction of the bind¬ 
ing energy, the kinetic energy of the initially free elec¬ 
tron increases in order to still fulfill the NEEC resonance 
condition. Since is inversely proportional to the 

squared electron momentum, this leads to a decrease of 
the NEEC cross section. 

Moreover, by comparing the cross section with the 
electron flux it can be seen that for Te = 350 eV and 
Te = 500 eV the resonance energies for M, N and O 
shell recombinations are located at the tail of the elec¬ 
tron distribution which shrinks exponentially. This be¬ 
havior is once more explicitly shown in Table IIIII where 
the resonant electron flux for Te = 350 eV and the in- 
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FIG. 2. (color online) (a) NEEC resonance cross sections Cneec 
for captures into the L, M, N and O shell (left axis) together 
with the electronic energy distribution (right axis), (b) Par¬ 
tial NEEC rate A^eec (l^ft axis) together with the correspond¬ 
ing probabilities of the available CSD (right axis). Results 
are presented for two plasma temperatures, 350 eV and 500 
eV— . Reproduced with permission from Phys. Rev. Lett. 
112, 082501 (2014). Copyright 2014 American Physical Soci¬ 
ety. 


main ingredients are the electron energy distribution in 
the plasma, the availability of capture charge states and 
the magnitude of the NEEC cross sections Cneec- In order 
to maximize the NEEC excitation rate one would prefer 
to have as many resonant free electrons as possible. As 
shown in Fig. [21(a) this is the case for the resonance en¬ 
ergies of NEEC into the L shell which lie in the region 
around the maximum of (j)e in the case of = 350 eV 
and Te = 500 eV. However, as can be seen from the 
charge state distributions shown in Fig. EJb), the num¬ 
ber of available charge states in the plasma will be limited 
for a given electron temperature T^. For instance, in the 
case of Te = 350 eV and Tg = 500 eV only Mo charge 
states up to g = 29 and g = 32 are present in the plasma, 
respectively. Opening the NEEC resonance channels for 
the L shell requires the Mo ions to be at least in the 
charge state g = 33. Hence, we conclude that L-shell 
capture is not possible for the considered temperatures 
in Fig. m 

In addition to the molybdenum CSD, Fig. j^Kb) shows 
also the nuclear excitation probability in dependence of 
the charge state g obtained from Eq. (HZI via a summa¬ 
tion over all contributing ou. Again the temperatures 
Te = 350 eV and Te = 500 eV are considered. Our re¬ 
sults show that the nuclear excitation probability grows 
with g because higher charge states open NEEC capture 
channels more deeply bound to the nucleus. 

Finally, looking at the dependence of Aneec on the elec¬ 
tron temperature Te we can conclude that increasing 
the temperature leads to a higher NEEC excitation rate 
for mainly two reasons: (i) higher temperatures involve 
higher available charge states which renders the capture 
into deeply bound electron shells possible; (ii) higher 
temperatures lead to an increase of the total number of 
available free electrons enhancing the resonant electron 
flux. 


D. Initial plasma conditions 


tegrated cross section—the so-called resonance strength 
5'neec—are presented for increasingly less bound capture 
shells. We have considered the example of Mo ions with 
initially fully occupied K, L, Mi, M 2 and M 3 shells, cor¬ 
responding to charge state g = 24. Moreover, NEEC into 
(^ 3/2 and c? 5/2 orbitals has been considered. The results 
clearly show that both the NEEC resonance strength and 
the number of resonant electrons decrease by going to 
higher principal quantum numbers. This behavior can 
be recovered for all capture channels and initial charge 
states, as long as the NEEC resonance energy is much 
larger than the plasma temperature. Hence, a “cut off” 
level can be found for each charge state g from which on 
the NEEC excitation can be neglected. Numerical results 
for Aneec shown in the following always employ such cut 
offs. 

For the hnal NEEC reaction rate calculation, the three 


In contrast to optical or infrared laser light, x-rays are 
able to produce directly inner shell holes by photoioniza¬ 
tion and penetrate much further into the material lead¬ 
ing to a very uniform irradiation of the solid. Moreover, 
the solid-state target is heated very rapidly by the laser 
pulse with duration of about 100 fs. On this time scale 
the ionic motion is negligible resulting in a plasma near 
to solid-state density. This rapid, isochoric heating of 
the plasma (the volume is nearly unaffected during the 
formation process) enables us to consider uniform tem¬ 
perature and density distributions at short times after 
the laser pulse. 

The laser-induced creation of the plasma is mainly 
dominated by two processes: (i) photoionization, and, 
(ii) Auger decay. In the latter process an electron on a 
higher atomic shell decays to an inner shell hole with the 
simultaneous emission of another electron into the con¬ 
tinuum. The process of photoionization prefers the inter- 
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action with inner shell electrons provided that the photon 
energy exceeds the ionization potential of the electrons. 
In the case of niobium, for instance, the ionization edges 
of the K and Li shells lie at 19 keV and 2.7 keV— , re¬ 
spectively. Hence, considering the resonant 4.85 keV x- 
rays, at most L shell holes can be produced. The laser- 
produced inner shell holes are then subsequently refilled 
by either radiative or Auger decay. 

In order to get a rough estimate of the initial plasma 
conditions directly after its creation, we follow the proce¬ 
dure used in Ref4i. We neglect the hydrodynamic expan¬ 
sion and the radiative losses for the first 100 fs when the 
laser pulse is present. (We note here that an extension of 
FLYCHK, the SCFLY code^, was particularly designed 
to model the plasma generation phase. Unfortunately, 
SCFLY is not freely available.) Assuming an instanta¬ 
neous equilibration time, equations of state (EOS) can be 
used to estimate the initial plasma conditions. Thereby, 
the deposited internal energy per pulse can be ap¬ 
proximated by e(J/ 5 ) = Ipee,h{W/crn'^)K{crn'^ / g)Tpuise{s) 
where k represents the photoabsorption coefficient. Us¬ 
ing the EOS tablesS^ e(Te,po); the initial plasma tem¬ 
perature immediately after the creation process can be 
determined. Note that the ionic density is thereby as¬ 
sumed to remain initially at its solid-state density po a-nd 
at room temperature. 

By using the LCLS laser parameters presented in Ta¬ 
ble H we obtain an initial temperature of k. 350 eV 
for the Nb target. For the calculation we used a photoab¬ 
sorption coefficient of 551.6 cm^/g. Analogously, the ini¬ 
tial stage of the Fe plasma can be estimated. Since due 
to the higher photon energy, the photoabsorption coeffi¬ 
cient (63.0 cm^/g) is much smaller in this case, we ob¬ 
tain a colder electron temperature from start on, namely 
Tefi 75 eV. 


E. Atomic effects in plasma expansion 

In dense plasmas, atomic processes such as recombi¬ 
nation and ionization are expected to play an important 
role on the plasma dynamics. In the following we sketch 
an approximate way to include the effects of the atomic 
processes on the hydrodynamic plasma expansion model 
described in Sec. Ill Di bv using numerical results from the 
FLYCHK codeia. 

FLYCHK was developed for describing plasma states 
like warm dense matter, highly transient states of mat¬ 
ter or extremely hot and dense matter. It employs a 
schematic atomic structure in order to provide a fast and 
widely applicable plasma diagnosis tool. Each atomic 
level is represented only by its principal quantum num¬ 
ber n. The atomic energy levels are computed from ion¬ 
ization potentials where the effect of the electronic con¬ 
tinuum depression occurring in plasmas is taken into ac¬ 
count by employing the model of Stewart and Pyatt2£. 
The population kinetics model implemented in ELYCHK 
is based on rate equations including radiative and col- 


lisional transitions between bound states, photoioniza¬ 
tion, collisional ionization processes. Auger decay, elec¬ 
tron capture, radiative recombination and three-body re¬ 
combination. These rate equations are solved for a finite 
set of atomic levels which consists of ground states, sin¬ 
gle excited states (n < 10), autoionizing doubly excited 
states and inner shell excited states for all possible ionic 
stages. 

With the model briefly described above, FLYCHK is 
able to determine ionization and level population distri¬ 
butions of a plasma (for some given conditions such as a 
given electron temperature and density). It can be ap¬ 
plied for low-to-high Z ions under most conditions of lab¬ 
oratory plasmas, in either steady-state or time-dependent 
situations^. 

The plasma under consideration here is a cold dense 
plasma. Considering again the numerical example of 
isomer triggering of ®^™Mo, the initial density of nio¬ 
bium ions is considered to be solid density and the initial 
electron temperature is several hundred eV. The result 
from FLYCHK shows that the time required to reach 
the steady state (with regard to the atomic processes) 
of such a plasma is on the order of several hundred fem¬ 
toseconds. With the radius of the considered plasma of 
around 1.8 /rm (according to the initial condition that 
the laser spot radius is around 1.8 /rm), the characteris¬ 
tic time scale for the plasma expansion is on the order of 
10 ps. Thus, the time scale for reaching the steady state 
(with regard to the atomic processes) is much smaller 
than the time scale of expansion. Therefore, we can as¬ 
sume that the steady state with regard to the atomic 
processes establishes at each time instant during the ex¬ 
pansion. As a first approximation, we can include the 
effects of atomic processes to the hydrodynamic expan¬ 
sion model by estimating the charge state of each time 
instant using FLYCHK. 

According to Eq. (1^ and to the ion density dynamics 
given by 

ni(i?) = n-po 

we may derive the dependence of the steady-state average 
charge Zi{R) of the plasma on the plasma size R using 
FLYCHK. In the equation above, rii o is the initial ion 
number density. Fig. [3] shows the average charge state 
of the ions for the case of ni_o = 5.5 x 10^^ cm“^ and 
Te^o = 350 eV. The average charge state decreases with 
reducing temperature and ion density. 

With the average charge state shown in Fig. [3] we 
can solve Eq. (1251) in order to study the plasma expan¬ 
sion. The results are shown in Figs. H) The notation 
tiexpn = dR/dt has been introduced. For comparison, the 
solution of Eq. (1331) for the case where the charge state is 
fixed to Zi Ri 23.86 is also presented. We observe that the 
expansion velocity is smaller when the effects of atomic 
processes in plasmas are taken into account. This is be¬ 
cause the average charge state, originally Zi « 23.86, 
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R/Rq 


FIG. 3. (color online). Steady-state average charge of nio¬ 
bium ions obtained from FLYCHK as a function of R (points) 
and its interpolation (solid line). We consider riifi = 5.5 x 
10^^ cm-® and Te.o = 350 eV. 


decreases due to atomic processes with the plasma ex¬ 
pansion. Thus, at larger times, the discrepancy between 
the two curves in Figs. |4] increases. We note however that 
the discrepancy remains small, on the level of 10%. Due 
to the plot scale, the discrepancy for is less visible, 
although it is exactly following the relation and it 
reaches the level of 10% at later times. 

In order to ensure the self-consistency of our approx¬ 
imation (i.e., that the steady state with regard to the 
atomic processes establishes at each time instant dur¬ 
ing the expansion), a time-dependent FLYCHK calcula¬ 
tion has also been performed, using the time-dependent 
Te shown in Fig. SJc), the time-dependent ion density 
given by Eq. (1^^ and R{t) in Fig. HJa) as input parame¬ 
ters. The result shows that Zi from the time-dependent 
FLYCHK calculation only slightly deviates (~ 5%) from 
the results in Fig. [3] for large i?, while the results agree 
well for small R. We conclude that on the degree of ac¬ 
curacy required for our calculations, the approximation 
used performs sufficiently well. 

We further study the behavior of the average charge 
state and the plasma radius for several initial tempera¬ 
ture values for both niobium (Figs. [5]) and ®^Fe (Figs. [6]) 
plasmas. Higher initial temperatures lead to higher 
charge states, with a very similar qualitative behavior 
at later times in the expansion when the temperature 
cools down. 


F. Total NEEC excitation 

Due to the hydrodynamic expansion of the plasma, 
the electron temperature and density decrease with time 
during the plasma expansion leading to a time-dependent 
net NEEC rate. Fig. [7] presents the time dependence of 
Aneec for Several initial electron temperatures Te.o- The 
ion density is always assumed to be at its solid-state value 
immediately after the plasma creation. As expected, in 
all cases Aneec drops down to zero with time, i.e., with the 
plasma expansion. With increasing initial temperature, 
Aneec takes larger values. 

The left column of Fig. [7] presents the NEEC reaction 



Time [ps] 

FIG. 4. (color online). Plasma radius R (a), expansion ve¬ 
locity Uexpn (b), and electron temperature Te (c) as func¬ 
tions of the expansion time. We use the parameter set 
Tii^o = 5.5 X 10^^ cm”® (niobium ions) and Te.o = 350 eV. 
Dashed orange curve: the result for the case where the steady- 
state average charge Zi{R) provided by FLYGHK and shown 
in Fig. [3] is included in the calculations, i.e., the atomic pro¬ 
cesses effects are taken into account. Solid blue curve: the 
result for a fixed charge state Zi « 23.86. 

rates A„eec for the case of ^^'"Mo triggering. Compar¬ 
ing the orders of magnitude of A„eec it can be seen that 
the NEEC excitation is strongly dependent on the ini¬ 
tial plasma conditions dictated by the laser parameters. 
Going from = 200 eV to Te.o = 500 eV enhances 
the NEEC excitation by approximately 7 orders of mag¬ 
nitude. The integration of Aneec over time provides the 
excited state occupation number per nucleus (compara¬ 
ble to piro, 

nOO 

PT"= * Aneec (30) 

Jo 

Performing the integration for the case of ®®™Mo re¬ 
sults in p"f'= Ri 1.8 X 10-20 for = 200 eV, p"f‘= a: 
1.4 X 10-15 for Te,o = 350 eV and p"|®‘= a: 1.6 • IQ-i® for 
Tefi = 500 eV. For all the three cases, the integration 
converged after approximately 6 ps. This value gives us 
the dominant NEEC interaction time during the hydro- 
dynamic plasma expansion. We recall that in Refi^, the 
plasma conditions were considered constant for 100 ps. 
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FIG. 5. (color online). Steady-state average charge of nio¬ 
bium ions (points) obtained from FLYCHK (a), and plasma 
radius i? as a function of the expansion time (b), for initial 
temperatures Te.o = 200 eV, Te^ = 300 eV, Te,o = 350 eV, 
Tefi — 400 eV, and Te,o = 500 eV. Numerical interpolations 
are shown by the curves. We use = 5.5 x 10^^ cm“®. 
The lowest charge state points correspond to cooling down 
the plasma to 50 eV electron temperature. 



Time [ps] 


FIG. 6. (color online). Steady-state average charge of ®^Fe 
ions (points) obtained from FLYGHK (a), and plasma ra¬ 
dius i? as a function of the expansion time (b), for initial 
temperatures Te.o = 50 eV, Te.o = 75 eV, Te^o = 100 eV, 
Te^o = 150 eV, and Te.o = 200 eV. Numerical interpolations 
are shown by the curves. The value ni,o = 8.5 x 10^^ cm“® 
was used. The end points of the curves correspond to the 
electron temperature cooled down to 15 eV for — 50 eV, 
25 eV for Te.o = 75 eV and Te.o = 100 eV, and 50 eV for 
Te.o = 150 eV and Te.o = 200 eV, respectively. 
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FIG. 7. (color online). NEEC reaction rate as a function of 
time for (left graph column) and ®^Fe (right graph 

column). In the case of initial electron temperatures 

of 200 eV, 350 eV and 500 eV have been considered; for ®^Fe, 
the corresponding values are 75 eV, 100 eV and 200 eV, re¬ 
spectively. 


Our present results show that this is not the case and 
significant changes in the plasma affect the magnitude of 
the secondary nuclear excitation after already few ps. 

The signal photon rates corresponding to the NEEC- 
induced triggering of ^^"‘Mo can be obtained by replacing 
Pee with in Eq. (051) . Assuming LCLS laser parame¬ 
ters, the signal rates evaluate to « 5.5 x 10“^"*^ s“^ 
for Te.o = 200 eV, i?"®®® « 4.3 x lO"® s”! for Te.o = 
350 eV and i?"®®® a; 4.9 x 10"^ s"! for Te.o = 500 eV. In 
order to fully exploit the repetition rate of the XFEL, a 
tape-station system must be used. Moreover, instead of 
having a single niobium target, a stack of several foils can 
be employed. This is possible as long as the transmitted 
laser intensity behind each foil is still strong enough to 
allow for the plasma formation. In this way an increase 
in the signal rates of approximatively one order of magni¬ 
tude may be achievable. Nevertheless, the isomer density 
of 10^® cm“^ remains the main limitation for obtaining 
higher count rates. 

We turn now to the NEEC reaction rates for ^^Ee pre¬ 
sented in the right column of Fig. [T] While the relevant 
time scales are similar to the ones of ®^™Mo, on the or¬ 
der of 6 ps, the magnitude of the NEEC rates is in com¬ 
parison tens of orders of magnitude below, and for all 
practical purposes negligible. Let us investigate the rea¬ 
sons for such a dramatic difference. Iron has Z = 26 and 
correspondingly smaller binding energies than molybde- 
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num [Z = 42), while the nuclear excitation energy is 
larger, 14.4 keV vs. 4.8 keV. The typical capture orbitals 
for both types of ions he in the M atomic shell. Cor¬ 
respondingly, the continuum electron energies required 
for NEEC will be much larger for ^"^Fe than for ^^’"Mo. 
While the M-shell NEEC requires continuum electron en¬ 
ergies between 3 and 4 keV, where the electron flux still 
has large values, in the case of ®^Fe the resonant contin¬ 
uum electrons should have more than 13 keV. For this 
energy value, the plasma at 75 eV temperature provides 
virtually no electrons at all, leading to the infinitesimal 
NEEC rate values in Fig. [T] 


G. Comparison between direct and secondary excitation 

In order to obtain comparable results for and 

Pel®® we hx the initial plasma conditions according to the 
chosen laser parameters. In particular, the laser inten¬ 
sity plays here an important role as outlined in Sec. lIII Dl 
Restricting ourselves to the LCLS laser parameters, an 
initial electron temperature of Te^o = 350 eV is estimated 
for the niobium target doped with ^^'"Mo isomers. Com¬ 
paring the corresponding excited state occupation num¬ 
bers and pII®®, we conclude that the indirect NEEC 
process is still about five orders of magnitude more effi¬ 
cient than the direct photoexcitation. One important as¬ 
pect is here that the time of interaction is about 60 times 
longer for NEEC than for the direct photoexcitation. In 
contrast to ®^™Mo, the indirect NEEC process plays no 
role for the excitation of ®^Fe, since the NEEC reaction 
rate values are for all practical purposes negligible. 

One important remark here is that the direct photoex¬ 
citation can be switched off by tuning the laser off res¬ 
onance. In contrast, the secondary excitation remains 
present since the XFEL-produced plasma is not sensi¬ 
tive to small detunings. In the case of ^^'"Mo for in¬ 
stance, the nuclear transition energy is only known up 
to an uncertainty of 80 eV. Since the XFEL bandwidths 
are several eV, this uncertainty can be a limiting fac¬ 
tor for the experimental implementation of the direct 
photoexcitation channel. The secondary NEEC process 
occurs instead in a plasma environment with a broad 
electron distribution where many resonance channels can 
contribute. Hence, it is much more robust than pho¬ 
toexcitation against such uncertainties in the nuclear pa¬ 
rameters. We also note that since the electron fluxes 
in the plasma are at the studied plasma temperatures 
several orders of magnitude larger than the correspond¬ 
ing photon fluxes, we neglect in our discussion the pos¬ 
sible secondary nuclear photoexcitation process in the 
plasma. Other studies2i“i2^ on hot astrophysical plasmas 
show that this may change at higher, keV plasma tem¬ 
peratures. The study of NEEC in plasmas was so far 
restricted to astrophysical environments^ir— or optical- 
laser-generated plasmasS^ where no equivalent of the di¬ 
rect photoexcitation channel under investigation here ex¬ 
ists. 


The values for p||®® are strongly dependent on the ini¬ 
tial plasma conditions. Theoretically, the initial plasma 
conditions and in particular the electron temperature 
should be controllable by changing the laser intensity. 
Since p||®® grows exponentially with Te_o saturating for 
large Te^o in the considered temperature range and 
furthermore depends linearly on the laser intensity / be¬ 
tween 10^^ and 10^® W/cm^, we expect for the consid¬ 
ered parameter range an exponential dependence of the 
NEEC-induced excited state population p||®® on /. We 
note that for this estimate we consider the same ini¬ 
tial size for the plasma, i.e., a constant laser focal spot. 
The direct photoexcitation in contrast follows pllf®® oc I. 
In addition, the XFEL intensity and consequently the 
plasma temperature can highly fluctuate from shot to 
shot. It is therefore reasonable to look at initial tem¬ 
peratures varying around the estimated value as it is 
realized in Fig. [T] For ^^"‘Mo, the indirect process is 
still of the same order of magnitude as the direct one for 
Tefi = 200 eV. For an initial temperature of 500 eV the 
indirect NEEC process becomes even competitive with 
the photo-induced excitation via the coherent XFELO 
^laser ^ ^ ^ 10“^^ although for the latter a 10 times 

longer pulse duration than for the LCLS has been consid¬ 
ered. In the case of ®^Fe, the indirect excitation channel 
remains negligibly small in comparison to the direct pro¬ 
cess even by considering an initial electron temperature 
of 200 eV. 

The reason why the nuclear excitation of ^^'"Mo and 
®’’Fe behave inversely with respect to the direct and in¬ 
direct excitations is two-fold. First, in the case of ^^'"Mo 
triggering the microscopic NEEC cross section is much 
larger than the one for resonant photoexcitation. This 
situation is reversed in the case of ®^Fe. The second 
and most important reason is related to the plasma con¬ 
ditions. In the case of resonant ®®™Mo triggering, the 
photoabsorption of the 4.85 keV photons by the niobium 
target {k = 551.6 cm^/g) is much stronger than the anal- 
ogon for 14.4 keV photons in iron {k = 63.0 cm^/g). Nat¬ 
urally, the more photons are absorbed, the more energy is 
initially deposited in the target leading to a higher initial 
electron temperature. Then, for NEEC it is important 
that the available charge states in the plasma and the 
energy of the free electrons are not that far away from 
the resonance condition. While this is the case for ^^'"Mo 
with a 4.85 keV nuclear transition energy, it is especially 
hard to realize in the limit of high nuclear transition en¬ 
ergies like for ®^Fe. 

Let us generalize our results based on the arguments 
presented above. A hrst insight into the relation between 
the resonant photoexcitation and microscopic NEEC 
cross sections can be gathered from the corresponding 
IC coefficients which are defined as the ratio between the 
IC and radiative decay rates. For high values of a^c we 
expect NEEC to dominate over photoexcitation, which 
is typically the case for small nuclear transition energies 
in the keV region. Next, the prevailing initial free elec¬ 
tron conditions are strongly dependent on the occurring 




14 


x-ray atomic photoabsorption. The photoabsorption co¬ 
efficient grows with increasing Z of the material and de¬ 
creasing photon energy. Moreover, shell edges may lead 
to a stepwise enhancement of k. However, a high elec¬ 
tron temperature simultaneously leads to higher avail¬ 
able charge states which again require smaller kinetic 
electron energies in order to fulfill the NEEC resonance 
condition. These opposite trends need to be balanced in 
order to match the available electron energies with the 
NEEC resonance condition for the open capture chan¬ 
nels. Typically, in the temperature range of hundreds 
of eV the secondary excitation contributes significantly 
only for transition energies below few keV. 

More concretely, in Table llVl we present the list of low- 
lying nuclear transitions first collected in Ref.— with the 
corresponding transition energy, multipolarity and IC co¬ 
efficient Oic. Also provided are the photoabsorption coef¬ 
ficient /t, the deepest ionization shell, an approximation 
of the initial plasma temperature Te^Oi the most proba¬ 
ble charge state at this temperature Z\, and finally the 
most NEEC-advantageous available capture shell in the 
plasma. The latter two quantities have been calculated 
with ELYCHK. Eor orientation, also the approximate 
ionization potential of the most advantageous available 
capture shell is presented. Since ELYCHK calculations 
are limited to the range Z < 80, the entries Zi and the 
capture shell for |o^Hg and | 2 ^Pb have been estimated 
by taking Z = 79. The initial temperatures in the table, 
with the exception of ^^'"Mo and ^’’Ee, discussed in de¬ 
tail in Sec. IHIDl have been obtained by accounting for 
the laser energy deposited into the sample with the help 
of mass photoabsorption coefficients^^. Due to unavail¬ 
able EOS tables, the averaged electron temperatures have 
been then roughly estimated by further considering en¬ 
ergy conservation of the inner shell photoionization and 
the first sequence of Auger decays. 

Based on the present analysis, we conclude that iso¬ 
topes with high nuclear transition energies like 
g 2 ®Sm and 50 ® Sn, which in addition only allow for NEEC 
into the M or even higher shells, will not present signifi¬ 
cant secondary nuclear excitation in the XFEL-produced 
plasma. Nuclei with small transition energies as it is the 
case for the first table entries down to lyRb are very likely 
to present significant secondary nuclear excitation in the 
plasma, due to the abundant electron flux at the required 
energies and the encouraging IC coefficient. As for the 
intermediate region with 6 keV< <13 keV, here a 
more careful analysis is required. The capture into L- 
shell orbitals in the case of gg'^Cs and 32 Ge speaks for 
more available free electrons in the plasma at the reso¬ 
nance energy. However, due to the low Z of germanium, 
it is likely that only in the case of 55 ^Cs secondary NEEC 
plays any role in the net nuclear excitation. In the cases 
of gg®Tm, 7 g^Os, gg^Tm and gf^’Ea, the rather large tran¬ 
sition energies, together with M- and iV-shell capture 
orbitals and the small IC coefficients indicate unlikely 
strong NEEC influence. Finally, the small IC coefficient 
in the case of yf^Ta also rather speaks against strong 


secondary NEEC. However, since the arguments above 
do not relate to the time scale difference between NEEC 
and photoexcitation and the unavoidable detuning for 
the photoexcitation channel, we conclude that dedicated 
simulations are required to draw solid conclusions for the 
nuclei in the intermediate energy region. 


IV. SUMMARY AND CONCLUSIONS 

The resonant driving of nuclear transitions with the 
XFEL is most likely to employ high-density, solid-state 
nuclear targets in order to improve the excitation rates. 
However, in this case secondary nuclear excitation in the 
produced cold, dense plasma might occur. We have quan¬ 
tified here the magnitude of the secondary excitation via 
NEEC taking into account the plasma dynamics after the 
laser pulse. We find that while photoexcitation can only 
occur during the pulse duration, on the order of 100 fs, 
the plasma effects last for several ps until plasma expan¬ 
sion diminishes completely the NEEC rate. In order to 
describe the time-dependence of the plasma parameters 
we employ a hydrodynamical model for the plasma ex¬ 
pansion together with atomic processes input which de¬ 
termines the ion charge distribution at each time instant. 
Plasma parameters such as initial temperature, charge 
state distribution or electron flux play a crucial role for 
the magnitude of the occurring secondary excitation. 

Our results show that for small nuclear transition ener¬ 
gies and advantageous free electron energy distributions, 
as it is the case for isomer triggering of ®^™Mo, secondary 
NEEC may exceed by as much as five orders of magnitude 
the direct photoexcitation by the XFEL. Our more realis¬ 
tic plasma condition analysis thus confirms qualitatively 
the results in Ref.—. Nevertheless, the signal photon 
rates for the ®^™Mo activation appear to remain below 
experimental detection, despite several possibilities to in¬ 
crease the triggering efficiency. In addition to enhanc¬ 
ing the number of irradiated isomers, one could envisage 
choosing XFEL parameters such that the conditions for 
NEEC triggering are optimized. The frequency of the 
XFEL radiation, for instance, has been so far chosen res¬ 
onant to the nuclear transition under investigation, in or¬ 
der to allow for direct photoexcitation. However, by fully 
exploiting the theoretical framework derived in this work, 
an optimized parameter set for more efficient secondary 
NEEC-induced triggering of ^^'"Mo may be deduced. 

On the other hand, we also show that for larger nuclear 
transition energies, in the case of the 14.4 keV Mossbauer 
transition in ®^Fe, secondary NEEC can be safely ne¬ 
glected. We may conclude that nuclear quantum optics 
experiments with bulk isotopically enriched ^^Fe samples 
in normal incidence will not suffer from strong decoher¬ 
ence rates due to plasma-related processes. Caution is 
however advised for extending these conclusions for ex¬ 
periments with x-ray thin-film planar cavities containing 
®^Fe iron layers. In these setups which operate in graz¬ 
ing incidence, the working principle of the x-ray cavities 
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Nuclide 

En 

E/ML 

K 

ion. shell 

Te,0 

capture 

E- 

-^lon 


^\C 



[keV] 


[cmVg] 


[eV] 

shell 

[keV] 

El 

Ml 

E2 

201-rr„ 

80 ng 

1.565 

Ml 

1975 

Ni 

670 43 

Ni 

0.378 


3.5e-h01 


rfPt 

1.642 

Ml 

1628 

Ni 

510 40 

Ns 

0.519 


1.2e-h02 


ifPb 

2.329 

E2 

911 

Ni 

790 45 

Ns 

0.645 



1.8e-h07 

gl^Sm 

4.821 

Ml 

362 

Ml 

740 42 

Ms 

1.420 


6.2e-h00 



4.850 

E2 

552 

Li 

350 24 

Ma 

0.410 



2.5e-h04 

e^Tm 

5.036 

Ml 

442 

Ml 

850 46 

Ml 

1.515 


1.6e-h00 


3?Rb 

5.236 

Ml 

328 

Li 

1080 32 

Li 

2.065 


7.5e-h01 



6.237 

El 

300 

Ml 

520 36 

N2 

0.465 

7.7e-01 



69® Tm 

8.410 

Ml 

118 

Ml 

630 43 

Ms 

1.468 


1.7e-01 


rTOs 

9.756 

Ml 

107 

Ml 

510 37 

Ns 

0.468 


4.8e-01 


ei'Tm 

10.400 

Ml 

291 

Li 

220 25 

Ni 

0.180 


3.1e-02 


sfLa 

10.590 

Ml 

166 

Li 

580 37 

Ms 

1.123 


3.7e-01 


i^Cs 

11.244 

Ml 

129 

Li 

1420 43 

L3 

5.012 


1.2e-h00 


SGe 

13.285 

E2 

124 

K 

220 18 

L3 

1.217 



4.1e-h02 

iePe 

14.413 

Ml 

63 

K 

75 12 

Ml 

0.093 


l.Oe-01 


IfSm 

22.507 

Ml 

28 

Li 

350 28 

Ms 

1.080 


3.1e-03 



23.871 

Ml 

12 

Li 

120 14 

N 2 

0.089 


l.le-02 



TABLE IV. Low-lying nuclear transitions suitable for resonant photoexcitation via XFELs. Transition energies and dom¬ 
inating multipolarities EjML, photoabsorption coefficients k, deepest ionization subshells, estimates of the initial plasma 
temperature Te^o, the most probable charge states Zi in the plasma, the most advantageous capture orbitals for NEEC avail¬ 
able in the plasma together with the corresponding ionization energies Eion and IC coefficients aic, respectively, are presented 
in the columns. The list is sorted by increasing nuclear transition energies En- 


is particularly sensitive to radiation damage. Based on 
our present results, a general set of criteria for identifying 
the parameter regime for which secondary effects in the 
plasma are important is presented. We anticipate that 
our findings will be of relevance for the layout of first 
nuclear excitation experiments at XFEL facilities in the 
near future. 
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